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Abstract 

Let X be a Harris recurrent strong Markov process in continuous time with general 
Polish state space E, having invariant measure \i. In this paper we use the regeneration 
method to derive non asymptotic deviation bounds for 



P x (\J*f(X s )ds\>t^s] 



in the positive recurrent case, for nice functions / with /u(/) = (/ must be a charge). We 
generalize these bounds to the fully null-recurrent case in the moderate deviations regime. 
We obtain a Gaussian contentration bound for all functions / which are a charge. The 
rate of convergence is expressed in terms of the deterministic equivalent of the process. 
The main ingredient of the proof is Nummelin splitting in continuous time which allows 
to introduce regeneration times for the process on an enlarged state space. 

Key words : Harris recurrence, Nummelin splitting, continuous time Markov processes, special 
functions, additive functionals, deviation inequalities, deterministic equivalent for additive 
functionals. 
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1 Introduction 

Consider a Harris recurrent strong Markov process X = (X t ) t > in continuous time with in- 
variant measure /i, taking values in a Polish space E. 
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If the total mass of /x is finite, X is called positive recurrent, null-recurrent otherwise. In the 
case of positive recurrence it is well known that for certain functions / such that //(/) = we 
have a central limit theorem for 



(i.i) 



as n goes to oo, see for instance Hopfner and Locherbach [TJ] and Touati |28j . 

In the null-recurrent case, the re-normalization y/n of (II .ip has to be changed. For that sake we 
have to consider what is called deterministic equivalent of additive functionals. The determin- 
istic equivalent has been introduced for Markov chains by Chen [5] and then been generalized 
to the context of continuous time diffusion models by Loukianova and Loukianov [20J and to 
any continuous time recurrent Markov process by Locherbach and Loukianova in [16J. It is a 
deterministic function 1 1- > v(t) such that v(t) — > oo as t — > oo and such that for any integrable 
additive functional A t , 



for any initial measure n. v(t) can be defined as follows. Take any fixed positive special function 
g of the process having /j,(g) > (see definition 12.41 below for the exact definition of special 
functions, for strong Feller processes, any bounded function having compact support is special) 
and define 



Jo 

where t] is an arbitrary initial measure. Then the strong Chacon-Ornstein theorem implies that 
for any other special function g' and any other initial measure if, 



Hence the deterministic equivalent is unique up to a constant in the sense that for two choices of 
the deterministic equivalent, v and v', we have that lim^oo v(t)/V(t) = c, where c is a positive 
constant. In regular models, v(t) ~ t a l(t), where I is a function that varies slowly at infinity. 
For example, for Brownian motion in dimension one, we have a = 1/2. The explosion rate v t is 
in general slower than the ergodic (positive recurrent) rate t. 

In the null-recurrent regular case, we have convergence in law of 



toSo W a , where B is a one-dimensional Brownian motion and W a the Mittag-Leffler process 
of index a, < a < 1, for certain functions / having //(/) = 0, see Touati ([28]) and Hopfner 
and Locherbach ([H]) for similar results for martingale additive functionals. In the fully null- 
recurrent case (i.e. null-recurrent, but not regular) no convergence in law holds true, which is 
a consequence of the famous theorem of Darling-Kac. 

Continuing in the spirit of the limit theorems of (II. ip and (II. 2p . the aim of this paper is to study 
deviation inequalities for additive functionals J * f(X s )ds for a certain class of centered functions 



lim lim inf P n (l/M < A t /v(t) < M) = 1 




E vUo 9(X s )ds) /2(g) 



t-+°°E v ,{f*g>(X s )ds) ti(g'Y 




(1.2) 
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/ (the class of bounded special functions that will be introduced below). More precisely we 
study the deviation bounds for 

pJ\^f(X s )ds\>^ +r >x^ 

for r\ > in the positive recurrent case, and in the general null-recurrent case, including the 
non-regular case, of 

p*(\j f(x s )ds\>v? + \y 

Thus we are interested in deviation inequalities in the moderate deviations regime. We will work 
in the most general situation, our results hold for any recurrent Markov process taking values 
in any Polish state space, without restrictions on the quality of recurrence. As a counterpart of 
this general approach, we have to restrict attention to the class of test functions which are the 
bounded centered specials functions. In the case of strong Feller processes, that means that we 
consider centered functions being of compact support. 

Our approach is based on the regeneration method. Regeneration times allow to split the 
trajectory of the process into i.i.d. excursions. In the one-dimensional case, such regeneration 
times are usually introduced as successive visits to recurrent points. Our aim is, however, to 
work in the frame of a general state space. In this general frame, recurrent points usually 
do not exist, and we use Nummelin splitting in continuous time, as developed in Locherbach 
and Loukianova [16], in order to overcome this difficulty This technique allows to introduce a 
recurrent atom for the process on an enlarged state space and thus to mimic the idea of recurrent 
points in higher dimensions. Even if the technical details concerning Nummelin splitting are 
somehow cumbersome, the evident advantage of this method is that it is conceptually very easy 
(regeneration means that we are able to work with i.i.d. variables) and that it works in any 
dimension and for any state space. 

For positive recurrent processes X, we obtain the following deviation upper bound. For any 
bounded special function / with //(/) = 0, we have 

P n (\ 1^ f{X s )ds\ > t^x^j < Ci exp (-C 2 t 2ri {x 2 A x)) + R t {x), (1.3) 

where 7r is any initial measure, for any < r\ < |. Here, C\ and C 2 are explicit positive constants. 
Rt(x) is a remainder term which is explicitly known and which is of order 0(exp (—\rt(x V 1))). 
In the general null-recurrent case, including the non-regular case, we have 

P n (j jf f(X s )ds\ > v(t)^ +ri x] < dexp (-C 2 v{tf {x 2 Ax)) +R t (x), (1.4) 

R t being of order 0(exp y—yjv(t){x V 1)J). Note that (jl.3p and (jl.4p yield Gaussian concen- 
tration in x. This is due to the fact that we consider functions / which are bounded and special. 
For this class of functions, the contribution of the additive functional up to the first regener- 
ation time can be controlled uniformly in the the starting point of the process, and it is this 
uniform dependence that is the basis of the Gaussian behaviour. In particular, for uniformly 
contracting Markov processes, Gaussian concentration is obtained for any bounded function / 
that need not be of compact support. 
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( II. 3p and (jl.4p represent a first and important step towards the study of moderate deviations in 
the fully general case of any recurrent process. The problem of obtaining such non-asymptotic 
bounds is of major importance for many applications. Let us cite just some of them : model 
selection or other non asymptotic problems for statistics of Harris processes, particle approx- 
imations of Gibbs measures, .... Our work is strongly motivated by applications to statistics 
(see for example [12]), in particular model selection, and it is important for such applications 
to obtain bounds that are valid for fixed t such that the constants involved are as explicit as 
possible. This is what we try to achieve in this paper. Note also that our results hold in any 
dimension and under any starting measure P w . 

Our deviation inequalities are stated in the moderated deviations regime. Note that the mod- 
erate deviations for additive functionals have been widely studied firstly in the case of discrete 
time, i.e. for Markov chains. We refer the reader to the work of Guillin Q12J), Djellout and 
Guillin ([8]), Chen and Guillin (jBj) for a survey on the subject. In the case of time continuous 
observations, less results are known. In the ergodic situation, Guillin and Liptser ([IB]) have 
studied the moderate deviations of f* f(X s )ds where X is a multidimensional ergodic diffusion 
and / any centered function belonging to L 1 ^). They use techniques of the stochastic calculus 
which are well-adapted to the particular diffusion case. Douc, Fort and Guillin ([9]) seem to 
be the first who deal with the general case of bounded additive functionals in the continuous 
time situation, but under assumptions implying sub-exponential ergodicity of the process (and 
in particular, positive recurrence) and assuming irreducibility of some skeleton chain (i.e. as- 
suming minorization of P m for some m on a petite set). All these results have been achieved in 
the positive recurrent case, specifying the quality of ergodicity - whereas in the present paper 
we work in the fully general recurrent but not necessarily positive recurrent case. 

Finally, let us cite a recent work by Loukianova et al. [TS] dealing with the same type of questions 
in the framework of positive recurrent diffusions in dimension one. In this case regeneration 
times are hitting times, and a precise control of these hitting times is provided in [18J. These 
kinds of results are not available within the general frame of the present article. 

The paper is organized as follows. Since we are using heavily the method of Nummelin split- 
ting and the concept of deterministic equivalent, section 2 gives a review of all known results 
concerning this technique that will be needed in the sequel. Section 3 gives the main results, 
for both positive and null-recurrent cases (Theorem 13. ip . In the case of null- recurrent but 
regular models, i.e. in the case when vt ~ t a for some < a < 1, some finer control of the 
Laplace transform of the length of the regeneration time yields better results than in the general 
null-recurrent case (see Theorem 13. 3p . Sections 4 and 5 are devoted to the proofs. Here, we 
also state a technical proposition which is a sort of generalization of Kac's formula, but more 
cumbersome in the context of Nummelin splitting (Proposition 14. ip . Finally, section [6] gives 
an idea of possible applications of our results in the framework of some interacting particle 
systems. 



2 Notation 

Consider a probability space (Q, A, (P x )x)- Let X = (X t )t>o be a process defined on (Q, A, {P x )x) 
which is strong Markov, taking values in a locally compact Polish space (E,S), with cadlag 
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paths. X = x Pj.— almost surely. We write L for the generator and {Pt)t for the semi group of 
X and we suppose that X is recurrent in the sense of Harris, with invariant measure //, unique 
up to multiplication with a constant. This means that for any set A e £ such that fi(A) > 0, 
limsup^^ lj^Xt) = 1 almost surely. Moreover, we shall write (J^tjt for the filtration generated 
by the process. 

We impose the following condition on the transition semi-group (Pt)t of X : 

Assumption 2.1 There exists a sigma-finite positive measure A on (E,£) such that for every 
t > 0, Pt(x,dy) — pt(x,y)A(dy), where (t,x,y) \-^p t (x,y) is jointly measurable. 

Example 2.2 1. In general, it is difficult to check whether a given Markov process is re- 
current or not, and the most used criterion for recurrence is the existence of a so-called 
Lyapunov- function for the generator of the process, see for example Meyn and Tweedie 

m w, «. ' ' 

We say that V G dom(L) is a Lyapunov-function, ifV > 1 and if there exists a constant 
a > 0, a constant b and a closed petite set C such that for all x, 

LV{x) < -aV{x)+bl c {x). 

The existence of Lyapunov- functions implies exponential ergodicity. This concept can be 
extended to obtain slower rates of convergence, still in the positive recurrent case, see 
Douc, Fort and Guillin JMj. 

2. In the context of interacting particle systems, (see Galves et al. (JTT$ and section^ below), 
recurrence can be shown via arguments using the dual process of the system. 

2.1 On the deterministic equivalent of additive functionals 

In the general recurrent not necessarily positive recurrent case, rates of convergence of additive 
functionals are given by what is called deterministic equivalent of additive functionals. This 
object has been introduced for Markov chains by Chen |5J and then been generalized to the 
context of continuous time diffusion models by Loukianova and Loukianov [20] and to any 
continuous time recurrent Markov process by Locherbach and Loukianova in [16] . In the context 
of one dimensional diffusion models, similar ideas to the notion of deterministic equivalent have 
also been developed in Delattre et al. ([I]). We start by resuming the most relevant results 
of (pL6j) on the deterministic equivalent that will be needed in the sequel. We first recall the 
notion of additive functionals. 

Definition 2.3 An additive functional of the process X is a M + — valued, adapted process A = 
{At)t>o such that 

1. Almost surely, the process is non- decreasing, right- continuous, having Aq = 0. 

2. For any s,t > 0, A s+t — A t + A s o 6 t almost surely. Here, 9 denotes the shift operator. 
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Examples for additive functionals are A t = J Q f(X s )ds where / is a positive measurable func- 
tion. Such an additive functional is said to be integrable, if fi(f) < oo. The deterministic 
equivalent of any integrable additive functional is a deterministic function v t— > v(t) such that 
v(0) = 0, v(.) is non-decreasing and v(t) — > oo as t — > oo. It satisfies that for any integrable 
additive functional A t , A t /v(t) is bounded and bounded away from zero in probability. In order 
to define the deterministic equivalent, we have to recall the notion of a special function (see 
also [H], g]): 

Definition 2.4 A measurable function f : E — >■ M + is called special if for all bounded and 
positive measurable functions h such that fj,(h) > 0, the function 



Note that in the case of strong Feller processes having locally compact Polish state space, any 
bounded function having compact support is special. 

By [16J, any special function g of X with fi(g) > defines a version of the deterministic 
equivalent via 



Jo 

for any arbitrary initial measure it. v(t) is called deterministic equivalent due to the following 
result (corollary 2.19 of |16j). 

Theorem 2.5 For any additive functional A of the process having E^(Ai) e]0,oo[, for any 
initial measure tt, we have 



Remark 2.6 The deterministic equivalent is unique up to a constant : for two choices of the 
deterministic equivalent, v and v' , we have that lim^oo v(t)/v'(t) = c, where c is a positive 
constant. 

In the sequel, depending on the situation, we shall fix a suitable choice of v(t). In the positive 
recurrent case, evidently v(t) = t, up to multiplication with a constant. In order to avoid too 
cumbersome notation, we sometimes also write v t = v(t). 

2.2 On Nummelin splitting in continuous time 

The proof of the deviation inequality is based on a very simple idea : the use of regeneration 
times that allow to divide the trajectory of the process into (almost) i.i.d. excursions. In the one- 
dimensional case, regeneration times are introduced as successive visits to recurrent points. For 




is bounded. 




(2.5) 
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Harris recurrent Markov processes with general state space, points are in general not recurrent. 
That is why we have to use the Nummelin splitting in continuous time, as developed in [16], in 
order to introduce a recurrent atom for the process. An atom is a set that, roughly speaking, 
behaves as a point for the process. Once a recurrent atom exists, we can introduce regeneration 
times that split the trajectory of the process into (almost) i.i.d. excursions. 

We recall briefly the construction of Nummelin splitting in continuous time. 

Introduce a sequence (<r n ) n >i of i.i.d. exp(l)-waiting times, independent of the process X itself. 
Let To := 0, T n := o~\ + . . . + a n and X n := Xt„- Then the chain X = (X n ) n is recurrent in the 
sense of Harris and its one-step transition kernel U 1 {x,dy) := J °° e~*Pi(a;, dy)dt satisfies the 
minorization condition 

U\x,dy) >al c (x)u(dy), (2.6) 

where < a < 1, fi(C) > and v a probability measure equivalent to //(• fl C) (cf [26], JO], 
proposition 6.7). The set C can be chosen to be compact. Note that the above minorization 
holds always under the only condition of Harris recurrence. We do not impose any further 
condition, neither irreducibility of some skeleton chain nor aperiodicity of the process. 



Remark 2.7 In some cases, the measure of assumption 2.1 A satisfies A ~ fi, and the densities 
Pt{x, y) are explicitly known, for example in the case of k— dimensional diffusions, under suitable 
regularity assumptions. If one can specify some set C and some time interval [s, t] such that 

inf p u (x,y) >0, A(C)>0, 

(x,y)eCxC,ue[s,i\ 

(w.l.o.g. also A(C) < 1) then Ii2.6\) holds true with 



a: 



\e~ s - e _ *l A(C) inf pJx, y) A 1 and v = A(- n C)/A(C). 

Ax,y)&CxC,u£[s,t] 1 



In particular, for multi- dimensional diffusions satisfying Hormander's condition on some set 
T, the classical results of Kusuoka and Stroock Iffl^ allow us to conclude that any choice of a 
compact set C cf will be possible. 



Then it is possible to define on an extension of the original space (fl, A, (P x )) a Markov process 
Z = (Z t )t>o, taking values in E x [0, 1] x E such that the T n are jump times of the process 
and such that under P x , ((Z}) t , (T n ) n ) has the same distribution as ((X t ) t , (T n ) n ). Here are the 
details of this construction. 

First of all, define the following transition kernel Q((x, u), dy) from E x [0, 1] to E : 

v(dy) if (x, u) G C x [0, a] 

Q((x,u),dy)= { T ^(U 1 (x,dy)-au(dy)) if {x, u) e Cx)a, 1] . (2.7) 
U\x,dy) iix^C 



We now recall the construction of Z t = (Zj, Zf, Zf) taking values in E x [0, 1] x E as given 
in [16]. Write w 1 (x, x') := J °° e^ t p t (x,x')dt. Let Z\ = X = x. Choose Zq according to the 
uniform distribution U on [0, 1]. On {Zq = u}, choose Zq ~ Q((x,u),dx'). Then inductively in 
n > 0, on Z Tn = (x, u, x') : 
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1. Choose a new jump time a n+ i according to 

-t Pt(x,x') 



dt on I. 



it (x, X ) 

where we define 0/0 := a/oo := 1, for any a > 0, and put T n+ i := T n + a n+ i. 

2. On {cr n+ i = £}, put Zj, +s := u, Zj, , s := x' for all < s < t. 

3. For every s < t, choose 

i p s (x,y)p t -s(y,x') 

Zt - +s R^o — A(dy) - 

Choose Zj- +B := for some fixed point Xq E E on {p t (z,^ / ) = 0}. Moreover, given 
Zr n +s = y, on s + u < t, choose 



7 1 Pu(y,y)pts-u{y ,x) , 

Pt-s{y,x') 



Again, on {p t _ s (y,x') = 0}, choose ^r n+s+u = x . 

4. At the jump time T n+ i, choose := Z\ n = x'. Choose Z^ n+i independently of 

Z s ,s < T n+ i, according to the uniform law U. Finally, on \Z\ +1 = u'}, choose 
Z*~Q((x', U i),dx"). 



Note that by construction, given the initial value of Z at time T n , the evolution of the process 
Z 1 during [T n ,T n+ i[ does not depend on the chosen value of Z\ n . 

We will write P v for the measure related to X, under which X starts from the initial measure 
ir(dx), and P„. for the measure related to Z, under which Z starts from the initial measure 
ir(dx) <g) U(du) ® Q((x,u),dy). In the same spirit we denote E n the expectation with respect 
to P n and the expectation with respect to F n . Moreover, we shall write IF for the filtration 
generated by Z, G for the filtration generated by the first two coordinates Z x and Z 2 of the 
process, and W x for the sub-filtration generated by X interpreted as first coordinate of Z. 

Write 

A := C x [0,a] x E. 
A is the recurrent atom of the process. Now we put 

So := 0, R := 0, S n+ \ := inf {T m > R n : Z Tm € A}, R n+ i := inf {T m : T m > S n+ i}. 

Then the sequence of IF— stopping times R n generalizes the notion of life-cycle decomposition 
in the following sense. 



Proposition 2.8 [Proposition 2.13 of ITh^ l 

a) Zft n+ . is independent of Ts n - for all n > 1. 

b) Zji n ~ u(dx)U(du)Q((x, u), dx') for all n > 1. 

c) The sequence of (Z Rn ) n >i is i.i.d. 
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Proposition 2.9 [Proposition 2.20 of IT5f] 

Let At be any integrable additive functional of X. Then, up to multiplication by a constant, for 
any initial measure 7r and any n > 1, 

E n (A Rn+1 - A Rn ) = MM) = EMi)- 
From now on, we shall fix a version \x of the invariant measure such that always 

/2(f) = E n / f(X s )ds. (2.8) 

JRi 

Moreover, we have the following: 

Proposition 2.10 [Proposition 4-4 °f WH] 

Let f be a measurable \l— integrable function. Put 



£„:= I*" f(X s )ds,n>l. 



Then the sequence (^ n ) n is a stationary ergodic sequence under P„. Moreover, for n > 2, ^ n is 
independent of Tn n _^. 

Remark 2.11 In the usual one- dimensional case, regeneration times R n allow to split in such 
a way that £ n is independent o/^ n _i. Our situation, however, is more complicated, since by 
construction, at a regeneration time, Z Rn depends on the state of the process one jump time 
before, i.e. on Zs n . This is due to the structure of continuous time and due to the use of Markov 
bridges (step 3. of the construction). 

Finally, let us recall the following useful result. 

Proposition 2.12 [Proposition 2.16 of ITty ] 
Let f be a special function of the process X. Then 

rRi 

C(f) :=su V E x / \f\(X s )ds < +oo. 

3 The deviation inequality 

In this section, we state our main result which is a deviation inequality for J * f(X s )ds, where 
/ is a special function of the process having = 0. Write 

N t := sup{n : R n < t}. 

In the following, we shall also use the following version of the deterministic equivalent. 

v* t = E u (N t ) + 1, (3.9) 

where v = C(Z Rn ) is given in (12.61) . The following is our main result. 
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Theorem 3.1 Let f be a bounded special function such that fi(f) = 0. Recall that C(f) = 
sup x Ex J 1 \ f\(X s )ds and put 

K(f) := ||/|U + C(f), B(f) :=max(K 2 (f),K(f)). 

1. Suppose that X is positive recurrent. Let 

F(X) = E(e- A(R2 - fll) ) 

be the Laplace transform of the length of a life cycle and put m = E(i?2 — Ri). Then for any 
u < m, 

A*(u) = sup [-Am - logF(A)] > 0, (3.10) 

A>0 L 

and we have the following results. 



(i) For any < r\ < |, for all x, for any initial measure it and for all t > 4m, 

p,(l/V,)*l>rf + '(f) 5+ V 

+8exp ( -t(x V 1)— A*(-m) ) . (3.11) 



4m 3 

(wj For a// a; and for any initial measure it, for all t > 4m, 

1 



P. (ijf /<*.)*! > yfix^j < 4exp (-^ (x 2 A x) 

+4e6XP (-6Jffc^) + 86XP {~ t(X V ^ A-( | m) ) ' (3 ' 12) 

2. Suppose that X is null-recurrent. Then there exists to, such that for all t > to, for all x, for 
any < r\ < |, /or any initial measure tt and for v\ as in Ii3.9\) . 



P n ( | / f(X s )ds\ > (v* t )^x ) < 4exp ( -^yy («?) V A x) 



+4eexp ( (u*)*^ + 8exp (-^KT(* VI)). (3.13) 



to defined through t£ o = 1. 



The bounds obtained in the above theorem are in general not optimal. But they are valid in a 
non-asymptotic framework. The quantities which are involved are almost all known or known 
up to some constants. Note also that as a function of the deviation level x, the above bounds 
exp(— Ct 2v (x 2 A x)) and exp(— C(v^) v (x 2 A x)) are of Gaussian type. This is due to the fact 
that we consider special functions allowing for exponential moments even in cases where the 
underlying process is not exponentially ergodic, see also Proposition 14.11 below. 



10 



Remark 3.2 In \3. y\) we do not use the usual form of the deterministic equivalent as chosen 
in \2. 5\) . But we have the following comparison result : 



1. For the choice of v t = E v f Q g{X s )ds as in h2.5\) . we have 

v t <C{g) + v(g)-v* t . (3.14) 
This can be seen as follows. Since g > 0, we have clearly that 



v, <^E,[l {Rn < i} I ^ L g{X s )ds\ 

n>0 ^ ' 
f rRn+i 

<C(g) + Y^^[kR n -i<t} / 9(X.)ds 

n>l ^ ^ Rn 



Using Markov's property with respect to J r R n _ 1 in the last expression and \2. 8\) . we obtain 
that 

v t <C(g) + ^g)[E v (N t ) + l}, 

and thus, 

v t <C{g)+/i{g)-v*. 
2. Moreover, we have the following lower bound. 

v t >fi(g)-v* t -2C(g). (3.15) 

Indeed, 

"Rn+2 \ \ rRN t +2 



v t > (l>,(l{*„<*} / ^ g(X s ))\ 

\n>0 V jR "+l ' JB 



g(x s )) 

R Nf 



> \E v (N t ) + l]n(g) - 2C(g) = v* t ■ fi(g) - 2C(g), 



since E„(f^ t+2 g(Xs))<2C(g). 



The above Theorem 13.11 can be improved in the null-recurrent but regular case. Here, in 
accordance with Theorem 3.15 of Hopfner and Locherbach, [H], we call a process regular, if for 
< a < 1 and a function I varying slowly at oo, the following is true. For any measurable and 
positive function g with < fi(g) < oo, we have regular variation of resolvants 

R 1/t g(x) = E x (J e-J s g{X s )ds\ ~ ^^(g) as t -> oo, (3.16) 

for fi— almost all x. Here we do not consider the case a = 1. We will show in proposition 15.31 
below that (13. 16[) is equivalent to the following : 

W(R 2 -Ri>x)~ % l ^ X \ as x ->• oo. (3.17) 
T(l — a) 
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Thus we are in the situation where R2 — R\ belongs to the domain of attraction of a stable law. 
In this finer control of N t is possible and it is well-known that 

v(t) ~ t a —— as t — > 00, 

see for instance Hopfner and Locherbach [TJ], theorem 5. 6. A. Moreover, 

1 _ p(X) ~ A Q /(1/A) as A -»■ 0, (3.18) 
see for instance Bingham et al. Q2], corollary 8.1.7). 
Then we get the following version of the deviation theorem. 



Theorem 3.3 Suppose that ( 3.11 ) holds. Let f be a bounded special function such that //(/) = 
0. Then there exists t > which is given explicitly in ( |5.^5| ) below, a function L varying slowly 
at infinity, such that for all x and for any initial measure it, for all t > to, for all < 77 < f , 

-Ptt 



jf /PQ^I > f3-*fc) < 4exp (-42^(7) ^ /(2_a) (^ 2 A x)L(t) 

^^ytf + ^ + 8exp (^-I^/( 2 -")(x V 1) A,-) 



+4e exp 



where 



Al = sup 



A>0 



\ogF(\)t a - a ^ ' 



L(t) 



At 1 " 5= 



> 



is positive and does not depend on t asymptotically : 

lim inf A* > (1 - a)a a/{l - a) > 0. 



t— >oo 



Remark 3.4 Note that in the above theorem, for rj = a/2, we obtain a rate of decay of the 
order of exp 

(_ t a/(2-a)y p ora = 2, this gives the rate of convergence exp(— t 1 / 3 ) which is better 
than the rate exp(— 1 1//4 ) obtained in Theorem \3.1\ item 2. (for rj = \). 

The above Theorems 13.11 and 13.31 do not use the Legendre transform of £ n = J^" i f(X s )ds 
since in general the law of the abstract regeneration times R n is not explicitly known. That is 
why the above deviation inequalities involve the constant K(f). In the following, we show how 
to compare K(f) and //(/). 



Remark 3.5 



1. By proposition 2.16 of ITb^ we know the following. Recall the definition of the set C of 
\2. 6}) . see also remark \2~1\ Let S be the first jump time of a Poisson process having rate 
lc(X s ), and let K be a constant such that for all x, 

E x I \f\{X s )ds < K, \f(x)\<K. 
Jo 

Then we have that 

IK 

C(f) <K + and hence K(f) <2K + . 

a a 
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2. We are now going to explore the relationship between C(f) and the invariant measure 
/J.{\f\) in some special cases. 

Suppose that the process X is strong Feller and positive recurrent. Let F be a compact 
set such that the support of f is contained in F. Let 

T c := inf{* >0:X t eC} 

be the entrance time in the set C. Recall that the measure v is concentrated on C. Then 



C(f) = wv% s / \f\{X s )ds< 

x&F Jo 



oo sup E X T C 

x&F 



Ri 



+ sup E x / \f\(X s )ds. (3.19) 
xeFnc Jo 



By continuity of the map x h-> E x J Q Rl \f\(X s )ds, there exists Xq G C fl F such that 



rRi 

sup E x / \f\(X s )ds = K X0 
xeFnc Jo 



Ri 



\f\(X s )ds. 



Moreover, there exists e > 0, such that for all x G B £ (xq), 



E x \ \f\{X s )ds> l -E XQ 



Hi 



\f\{X s )ds. 



Then we have that 



fi(\f\) = E u j Rl \f\(X s )ds > [ u(dx)E x [ VlpQtfo 

Jb £ (x ) Jo 

>~.v(B e (x ))-E X0 J** \f\(X s )ds. 
Putting together Ii3.19\) and Ii3.20\) . we conclude that 



C(f)< 



oo sup E X T C 

x£F 



A\f\)- 



(3.20) 



(3.21) 



HBe(xo))' 

3. Suppose that X is a positive recurrent one-dimensional diffusion process 

dX t = b(X t )dt + a(X t )dW t . 

In the case of dimension one, we can avoid Nummelin splitting since successive visits of 
recurrent points allow to split the trajectory into i.i.d. excursions. More precisely, let 
a < b be two recurrent points of X and define a sequence of stopping times (S n ) n , (R n ) n 
as follows. Sq = Ro = 0, 



Si = inf{* > : X t = b}, R x = inf{t > 5i : X t = a}, 



(3.22) 



and for any n > 1, S n+ i = R n + S\ o 9n n , R n +i = R n + Ri ° 6n n - Using this sequence 
(Rn)n, under additional regularity assumptions (see /7^/ for the details), we have that 

C(f) < Ml/1), 

for any function f having compact support. Here, the constant k is explicitly known. 
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4 Proof of Theorem 



The following proposition will be crucial in the sequel. It states that additive functionals built 
from bounded special functions admit a certain number of exponential moments. 



Proposition 4.1 Let f be a bounded special function. Recall that K(f) 
we have for any initial measure ir and any n > 1, that 

In particular, we obtain for any < A < K(f)^ 1 , 

E.(e A ^) < i + XPK UY = Y^y h 

P>1 

and if fi(f) = 0, 

^) < i + g \PK(fy = i + f K ^ f) ■ 



E, 



, + C(f). Then 
(4.23) 



Proof Evidently, we have that 

"Rn fRn 



e n = / •■•/ f(x tl )...f(x tp )dt 1 ...dt p 

J Rn-l J Rn — 1 



< 



"Rn l-Rn 

P! / ••• / l{ti<...<^}|/|(^ti)...|/|(^O d *l--- d V 

I Rn-l JRn-l 



Taking expectation and conditional expectation with respect to Tt v -n we get 

l-Rn fRn 

K ■■■ l {tl <...< tp} \f\(X tl )...\f\(X tp )dt l ...dt p 

J Rn-l J Rn-l 

l-Rn fRn 

= K ■■■ hu<...<t p . l} \f\(X tl )...\f\(X tp _ 1 ) 

J Rn-l J Rn-l 

Ez. . I \f\(X s )ds 



z * P -i 



o 



dt\ . . . dtp_\. 



But note that for any fixed z, for M 



\f\(X s )ds 



rRi 

< ME 2 (T\)+E 2 / \f\(X s )ds 

JTi 

< ME z (Ti) + E z (E z i [ l \f\(X s )ds) 

1 Jo 

< ME^+Cif). 

(Compare to proposition 2.16 of [IS]-) But by construction, for z = (x,u,x'), 

ttji (rp \ f°° , -tPt{x,x') 

E 2 (T\) = / te t — -dt, 
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(4.24) 



(4.25) 



and this expression does only depend on x and x' . According to proposition 4.1 of [T7], we have 
that 

L[Z\\Z\ = x)(dx') = v}(x, x')A(dx'). 

Taking now conditional expectation in (I4.25P with respect to T tv _ 2 and putting all these results 
together, we obtain 

"Rn 



E. 
= E. 
< E 
= E. 



Rn 
Rn-1 



Rn- 



R n 



Rn—l J Rn- 

Rn fRn 



l {tl <...< tp} \f\(X tl )...\f\(X tp )dt 1 ...dt p 
l { i 1 <...< w} |/|(X tl )...|/|(X^ 1 )E Zfp _ i 



Ri 



\f\(Xs)ds 



Rn- 



Rn-1 J R 



Rn- 
Rn 



l {tl <...< tp _ l} \f\(X tl )...\f\(X tp _ 1 ) ME^JTJ + Cif) dt x ...dt v _ x 



|/|(X.)(ME z .(T 1 ) + C(/))ds 



Rn 



R 



n—1 J Rn—l 



Rn 



l{ti<...<tp- 2 }\f\{Xti) ■ ■ ■ \f\( x t p - 2 ) 

- rRi 

L/o 

l{t 1 <...<V 2 }l/l(^ 1 )---l/l(^- 2 ) 

Ri 



dt\ . . . dtp— 2 



E 



z t 



p-2 



\f\(X s )(M + C(f))ds 



where the last equality follows from 

-Ri 



dt\ . . . dt p ^2i 



E z ^ * IfKXjEz^Tjds = E z J™ l {s<Rl} \f\(X s ) Qu\X,,x?)A(dJ)E {Xtt * ) ( k T 1 ?jds 

POO 

= E z / l {s<Rl} \f\(X s )ds, 
Jo 



since 



u\x,x')Mdx') te- t Pt ^ ,X '\ dt= 1. 
1 ^ 1 'Jo u\x,x<) 

Taking successively conditional expectations with respect to J-j 3 , ■ ■ ■ , J^ti yields the result. 
A first step in order to prove Theorem 13.11 is the following proposition. Let 



v t 



—t if X is positive recurrent 
vl if X is null-recurrent 



Proposition 4.2 Let f be a bounded special function such that fi(f) = 0. Then there exists t Q , 
such that for all x > 0, for any initial measure 7T, for any < 77 < |, < 5 < 2t7, for any fixed 
choice of a deterministic equivalent v t and for all t > to, 



f{X s )ds\>vt x) < 4exp(-- max{K2{f) jK{f)) v t 



x 2 A x 



2r/-S 



+4e exp 

Here, to is given by the equation v to = 1. 



' v} + \ ) + 4K(N t > v] +5 {x V 1)). 



QK(f) 



15 



Proof Write £ n := J R f(X s )ds } then the £„,n > 2, are identically distributed random 
variables having mean zero, such that £ n and £ n+ 2 are independent. It can be proven that for 
any p > 1, E(|£ n | p ) < oo, by the properties of a special function. More precisely, we have that 
for any n > 1, 

E(\U\ p )<P ] -K(f) p , (4.26) 

which will be shown in proposition 14. II below. Thus for A > 0, sufficiently small (A < K(f)^ 1 
suffices) , 

Z{\) := E^(e A? "), n > 2, (4.27) 
exists and is finite (and does not depend on n). 

Since the £ n are not independent, but only 2— independent, we have to proceed in the following 
way. Firstly, define a sequence ' by 

^> = {o" ei n eif° dd }- < 4 - 28 > 

Then define a second sequence ^ by 

t(2) _ f £n if n even 1 
Sn \ elseif J ' 



Note that £ n is not independent of J 7 r„_ 1 , but it is independent of TR n _ 2 . That is why we 
introduce the following two sub-filtrations, associated to the sum of odd and the sum of even 
terms. Let 

S« := cr{R k , : h <n,k odd }, 

and 

£i 2) := f < 2) :k<n } k even }. 

Moreover, let 



iV^ := sup{n : n odd , i? n < t}, iV 4 ( ' z '' := sup{n : n even , < t). 



(2) ._ 



Then it is immediate that + 2 is a (Gn') n — stopping time and N^ 1 + 2 a (Gn') n — stopping 
time. Moreover, we evidently have that < N t , < N t . 



(2) 



,(2)v 



Now, for A sufficiently small, let 



Ml := exp(A^^ 1) ) ■ Z{X)-^^\ M n 2 := exp(A^f) ■ Z(Xp 



(2)> 



[n/2] 



fc=2 



fc=2 



where [.] denotes the integer part of a real number. Then (M*) n and (M 2 ) n are discrete 
martingalges (^i 2 ^— martingales, respectively). Hence using Doob's stopping rule for pos- 
itive super-martingales we get 



exp(A £ ^ 1) )exp(-[(iVi 1) + l)/2]logZ(A)) 



(i) 



fc=2 



< 1 



(4.29) 
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and also 



exp(A d 2) )exp(-[(iVf ) + 2)/2]logZ(A)) 



k=2 



< 1. 



Now, we proceed as follows. Evidently, 

"-Rjv t +2 



where I = > v} +v x/3) 



7V, (1) +2 



/(X,)ds>0)</o + /i + J 2) 



7V t (2) +2 



A = **( E tf } ^ ^ + V3), / 2 = P.( E d 2) > ^ + ^/3). 



fc=2 



fc=2 



We start with a study of the first term J . Let 

Z(X) = E^ 1 . 



Then, for < A < K(f) 



-i 



\[6 > v} + \/Z] < exp 



•(^ + "As/3-logZ(A)) 



But due to (14.261) . we have that logZ(A) < ^^Wf) ' anc ^ thus, taking A = \K(f) 



h = > vr v xK(f)/3] < e exp 



1 A+v 



vr x 



Let us now turn to the study of ii and li- 
Note that for any k > 0, 



AT (1) +2 



h < K( ^ > v! +V x/3; N t <k)+ ¥ n (N t > k). 



k=2 



By g2SD, we have that for any < A < K(f)~ 1 , 



£(A) < i + 5>n/T- 



n>2 



Hence 



1-AC(/) 



l-AiC(/) 
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Hence for any < A < K(f) l , using (I4.29P and recalling that < N t 



N^ 1} +2 



7V t (1) +2 



K[ E & * v! + \/S ;N t <k}<F n [J2 > v! + \/3 ; iV« < k] 

k=2 k=2 

< ¥ w IM 1 



.^ (1)+2 > exp(A< + V3-log^(A)[(iV t (1) + l)/2]) ; iV t (1) < k 



< P. 



M^x, > exp(A^ + \/3 - /^l ^ + l)/2]) ; ivf } < fc 



M 



1 - XK(f) 
X 2 K(ff 



[{k + l)/2}) 



<»p(-(^Axi-^L*)). 



In the same way we get that 



Now, take 

Then we have that 



k = v] +s (s VI). 



p w [ E *£° ^ J ^ ^ v ^ ^ ex p *))) 



k=2 



(4.32) 



(4.33) 



7V t (2) +2 



P.[ E ^ - ^x/3;N t < v] +5 } < exp (-^ 1+<5 /i(t, x))) 



fc=2 



where 



/i(t, a?) := (x V 1) sup 



AQeAI) _ X 2 K(f) 2 \ 

l+S-ri 1 - \K(f) J 



(4.34) 



0<A<A'(/)-i ^3^2 

It can be shown, see for example Birge and Massart ([2]), lemma 8, pages 366 and 367, that 



sup Xy - 



2„,2 



X 2 v 



1 - Xv 



> 



0<X<l/v 

This is seen as follows. A simple calculus shows that 

AV 



2vy + 4t> 2 



(4.35) 



sup ( Xy — 

0<\<l/v 



1 - Xv 



(A*)V 
1 - AV 



X*y 



where 



A* = - ( 1 

v 



y + v 



< 
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It follows that 



l-X*v \\y + v ) yv + 2v 2 + 2v 2 (l + $f)V2 ' 
and using (1 + y) 1 ^ 2 < 1 + y/2, one gets the desired inequality. 
Using f04jl and p5| . we get for any i > i such that (^ )^ +5 ~' 7 > 1, 

Ht ' X) -3 (XV1) TO r <i^> + 12^(/) ^ « <436) 

Note that the right hand side of (I4.36P tends to zero at speed v^ 1+2ri ~ 25 . However, this right 
hand side has to be multiplied with vl +s , compare to f)4.33p . which yields the term vf 1 ^ 5 which 
does not tend to zero since by assumption, 5 < 2r\. Thus, together with (I4.3ip . 



Pirl/ f(X s )ds>v t x] < 2ex P 

Jo V 42 max(K 2 (f),K(f)) J 

1 4 



+e exp f -^_^\J + 2F n (N t > v] + \x V 1)). 

Applying the same argument to — / instead of / yields 

M |/ WW > -,*) < 4e*p (-- ^ (4.37) 







+2e exp ( --^L^ v ? + \ ) + 4P 7r (iV t > V 1)). 



Moreover, note that 



r--Rjv t +2 i \ / rRN t +i 1 i 

P ff (|/ /(X s )ds| > < pJ| / /(i s )^|>-o 



+Pj I / f(X s )ds\ > -vf\ 



and, using Markov's property with respect to JF t and to J-r n +1 , 

F n (\ j* N±+1 f{X s )ds\ > ^vf^xj < E T E Xt e x ^ lmx ° )ds e- v t +Vxx/2 . (4.38) 



Write 

-R 1 



Z(X)=E 7T E Xt (expX / \f\(X s )ds). 
Jo 

As before we have that 

~ XK(f) 
IogZ(A) < 



l-XK(f) 
19 



and thus, taking A = \K(f) 1 



LP, | | / f(X s )ds\ > -vf'x ) < e exp ( --^ v f" x 



In the same way we get that 

P tt | I 1^ f{X s )ds\ > ^ +v x\ < e exp {^^J^ v t +V x 



Remark 4.3 The most delicate point in the above proof is ( |^ . 36] ) which gives a lower bound 
converging to zero. It is important to be able control the speed of convergence of this expression 
to zero. If we decided to work with the Legendre transform, then the above proof remains true 
by replacing h4-34\ ) by 

(x A 1 
A-r- logZ(A) I . 
v t 2 V 

In this form it is evident that h(t,x) — > as t — > oo, and developing Z{\) = 1 + ^A 2 + o(A 2 
where 

-R2 



a 2 := E( / f(X s )ds)\ 
Jr.! 

would yield the same speed of convergence. However, using this approach we would not be able 
to put hands on the exact form of the constants appearing in the remainder term o(A 2 ), and 
that is why we decided to use the upper bound of log Z(X) as proposed in the proof. 



In order to make use of Proposition I4.2[ we have to control the number of life cycles before 
time t, N t . We first study the positive recurrent case. Recall that v t = —t in this case. 



Proposition 4.4 Suppose that X is positive recurrent. Then 

P 7r (iVi > v t (x V 1)) < 2 exp (--Vt(x V 1) suph t (X,x] 



2 



x 



where 



h t (X,x) := — (! — -7-77-7] logF(A) - A- "' 



t(xVl)J 2(x VI) 



and where 

F(X) = E(e~ x(R2 - Rl) ) 

is the Laplace transform of the length of a life cycle. Moreover we have for any u < m that 

A*(u) = sup [-Am - logF(A)j > 0. (4.39) 

A>0 L 
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Proof of Proposition 14.41 Writing k t := [v t (x V 1)], for any A > 0, 

F^Nt > v t {x V 1)) = F^{R kt <t)< F w (R kt - R 1 < t) < F^ X ^ R ^ > e- xt ). 

But, using the definition of £n and as in (14.281) . with / = 1, and the same technique as 
above, 

P^e-AG^-Ri) > e -At) < p^g-AE^d 15 > e -W) + P7r ( e -AE*t 2 d 2) > e -A*/2) 

< 2j p (A) (fc t -l)/2 e At/2 < 2j p (A) M^ eA4/2; (4 _ 4Q) 

where 

F(X) = E(e-^ R2 - R ^) 
is the Laplace transform of the length of a life cycle. So write 

ht{X,x) :=-(!- , 2 W1 , 1 logF(A) - A- ' 



v t (x VI)/ v y u t (a; V 1)' 
Then 

F^Nt >v t (xVl))< 2 e -W*vi)su PA MA,2)_ 

Let us finally show that 

A* (it) = sup \-Xu - logF(A)] > 

A>0 L 

for any u < m. This is evident using the fact that for A — > 0, — logF(A) ~ 1 — F(X) ~ mX. 
This concludes the proof. • 

Proof of Theorem 13. 1L item 1. The proof follows from Proposition 14.21 taking 5 = 
and from Proposition 14.41 in the following way : A straightforward calculus shows that, for all 
t > 4m, since x V 1 > 1, 

3/2 \ 3 2 

s\xph t (X,x) > -sup ( - logF(A) J = - A*(-m). 



This concludes the proof of item 1. of Theorem 13.11 • 

In the null-recurrent case, recall the choice of a deterministic equivalent v% of (13. 9p . Then the 
following proposition allows us to control the deviations for N t in the null-recurrent case. 



Proposition 4.5 



Proof 



E n (N!)< P \(v;r. 



N t = Y 1 { r «<*}' 

n>l 



and thus 



Nf=p\ <*}•■■■• 1 {R np <t} =P } - l {Rn p <t}- 

l<ni<n2<...<n p l<ni<n2<...<n p 
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Hence, using successively Markov's property with respect to J r n np _ 1 etc, we get 
K(Nf) < P \ E -( 1 {«>- 1 <*} E ^ > _ 1 (E !{**<*}))■ 

I<ni<n2<...<n p — i n p >l 

But note that for any n, 

fc>l fc>2 k>l 

= 1 + E U J2kR k <t} =v*(t), 

k>l 

since £(Z^J = z/. Induction on n then yields the assertion. • 
We get the following corollary. 

Corollary 4.6 We have that 

F n (N t > (v;) 1+s (x V 1)) < 2exp (~(vl) s (x V 1)^ . (4.41) 

Proof By proposition H~5*l 

E,(exp(AiV t )) < ' 



1-Xvf 

for A sufficiently small. Choosing A = {yl)~ l /2, we get that 

¥ n (N t > (v* t ) 1+5 (x V 1)) < e- A(l, ? )1+ ' (:cV1) E^exp(A^)) < 2e~^'^ 



Proof of Theorem 13.11 item 2. Choosing 5 = rj, item 2. follows immediately from Proposi- 
tion H2] and from (I4.4ip above. • 



5 Proof of Theorem 13.3 



Proof of Theorem 13.31 The proof is a slight modification of the proof of proposition 14.21 We 
go back to the proof of proposition 14.21 We fix the following choice of k (compare to (I4.32p ) : 

1 I- a 

k = k t = 2 + (x V l)t 7 7 . where 7 = a + 2rj- . (5.42) 

t(l/A t J I — a 

Here, \% — >• will be defined in (I5.44p below and I is the function of (I3.18p . 

As in (14 .40 p . we have for any A > 0, 

K(N t > k) < 2e -M- 1 °^W(fe-2)-A*] j 
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and therefore, due to the choice (I5.42p . we have to find 



sup 

A 



logF(A)0 VI)*" 



J(l/A t 



At 



Writing L(t) = 1(1/ Xf), this last expression can be rewritten as 



(x V 1) sup 

A 



1-7 

1 — I t 1 -" 



L(t) x V 1 



-'■I 



> t~(x\/ 1) A 



t ! 



where 



A* = sup 



logF(A)* Q ^i - Xt^ 



Lit) 



(Note that = 1 — ^-.) Write for simplicity 



2 n 



1— 



Then, 



Aj = sup 

A>0 



s t := *i- 



logF(A)-- S ?-A St 



* x "^ 



sup 

A>0 



2>7 



logF(A) 



/(1/A t 



and we have to show that A^ is positive. Note that since R2 — R\ does not possess any moments, 
we are not able to develop F(X) near in the usual way. But we can use (I3.18p . That's why 
we take A of the form A t — > at a speed that will be precised in (I5.44p below. Note that 

\ogF{X t ) = log [l - (1 - F{X t ))] . Due to ( Eg} we have that 



and hence 



-logF(A 
A*>-logF(A t 



W i(l/At) 
1 



s« - X t s t ~ A"s° - X t s t , 



l(l/X t 



a" — X t s t ~ A" - X t s t as * —7- 00. 



Maximizing A a — As t with respect to A suggests the choice 



X t = a 



1/(1-0,) -1 
5 t • 



(5.43) 



(5.44) 



For this choice we get that 



A" s t Q - X t s t = (1 - a)a a/(1 - Q) > 0. 



This implies that A* is strictly positive eventually and that 

lim inf A* > (1 - a)a a/{l ' a) > 0. 



t— >oo 



We continue the proof following the lines of the proof of proposition 14.21 (14.331) is true with 
v replaced by t% +v and v} +s by k t . Here, 



h(t, x) = sup ( A 



A>0 \ OK t 



3k f *• 



A 2 X(/) 2 \ / xAl X 2 K(ff \ 

-§-» i-xK{f)) 31 3M-f-" i-mf)) 
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where k t = ^ + P/L(t). Let 



t 7 

t such that for all t > t , J^ r ^~ v > L ( 5 - 45 ) 



Then for t > t , with B(f) = max(K 2 (f), K(f)), 

1 



h{t,x) > i (x 2 Al) L{tf 



[2(xAl)K(f) + 12K 2 (f)\ ~ 42B(f) 



£-2 7 +a+2r; > X 2 A 1 W^\2 ^-2-y+a+2ri 



Remark 5.1 Note that in the regular case, we even have convergence in law of 
(-i= J^f(X s )ds) t to aB o W a , where B is a one-dimensional Brownian motion and W a the 
Mittag-Leffter process of index a, i.e. the process inverse of the stable subordinator of index a, 
see for instance Touati 



There are various examples where the exact form of the Laplace transform is known. Brownian 
motion in dimension one is the most famous example. 

Example 5.2 Let X be the one- dimensional standard Brownian motion. In this case, we can 
define the regeneration times without Nummelin splitting in the following very simple way. Let 

R n = mf{t >S n :X t = 0}, S n = mf{t > R n ^ : X t = 1}, Rq = 0. 

Then _ 

F(X) = e- 2V2 ~\ 

see Revuz-Yor flWJjl). In this case it is possible to take v t — yi (see for example Hopfner and 
Locherbach (fT^j), theorem 5. 6. A.). 



We conclude this section with the following proposition showing that the notion of regularity 
(I3.17P is intrinsic of the process and does not depend on the concrete splitting we are using. 



Proposition 5.3 Suppose that for < a < 1 and a function I varying slowly at oo, the 
following is true. For any measurable and positive function g with < [i(g) < oo, we have 
regular variation of resolvants 



R 1/t g(x) = E x (J 



OO \ 1 



e-- s g(X s )ds ) ~ t a —fi(g) as t -> oo, (5.46) 



for fj,— almost all x. Then we have 



F(R 2 - Ri > x) ~ * l ( x \ as x ->■ oo. (5.47) 

r(i — a) 
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Proof By theorem 3.15 of [2], we know that ( I5.46P implies weak convergence 

{A tn ) t > ^^ i)r 



n a /l(n) 

for any additive functional A t of the process having < E^Ai) < oo. This convergence holds 
true in D(R + , R), under P n for any initial measure tt. W a is the Mittag-Leffler process of index 
a. 

From now on we fix v as initial measure. Interpreting X as first coordinate of Z and using 
Chacon-Ornstein's theorem, we have in particular weak convergence of (N tn )/{n a /l(n)), under 
¥ u . Now, write v(i) = t a /l(t) and let a(n) be its asymptotic inverse, i.e. a(n) ~ n 1 ' a l(n). Note 
that a(n)/n — >■ oo since a < 1. JVt is the inverse process of the sequence of regeneration times 
R n , which implies that 

R 

converges weakly as n — > oo, under P„, (5.48) 



an 



where the limit is non-degenerated. But R n = Ri + (R 2 — Ri) + . . . + (R n — R n -i)- We would like 
to deduce from this that necessarily (I5.47P is true. Unfortunately, due to the complex definition 
of the Nummelin splitting, the (Rk — Rk-i)k are not independent, but they have all the same 
law. Independence holds true only after "nearly" exponential times, the reason for this being 
the choice of the new jump times in step 1. of the splitting algorithm. 

That is why we proceed as follows. Define for any k > 0, Tf. = inf{T n : T n > R^} the first 
jump of the process after the k— th regeneration time. Then it is straightforward to show that 
under F u , Tk — Rk is exponentially distributed with parameter 1, and that the random times 
Tk — Rk, k > 0, are independent. Hence, 

^2_o(Tfc — . . Yll-niTk — Rk) n 

fe ~ u -> 1, which implies that fe ~ u -)■ 

n ' a(n) 

almost surely. This together with (I5.48[) yields weak convergence of 



a(n) 



(5.49) 



But Rk+i — Tk, k > 0, is an i.i.d. sequence of random variables under F u , thus (I5.49P implies 
that necessarily the law of Rk+i — Tk belongs to the domain of attraction of a stable law, see 
for instance Feller, [TO] . XIII. 6. This implies that 

F U (R 2 - fi > x) ~ x~ a L(x) 

for some function L varying slowly at infinity. Since 

¥ U (R 2 -f 1 >x)< ¥ U (R 2 -R 1 >x) < F U (R 2 - Ti > x/2) + P„(fi - Ri > x/2) 

= F U (R 2 - fx > x/2) + e~ x '\ 

we deduce that 

x \—¥ P^(i?2 — Ri > x) varies regularly at infinity with index a. 
This implies (15.4=7)) . • 
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6 Application to some interacting particle systems 



In some applications it is interesting to know at which speed empirical means | f(X s )ds 
converge to the - in general unknown - invariant measure n(f)- This is most often the case in 
statistical applications. 

Consider for instance the following interacting particle system which has been studied in Galves 
et al. [11]. Particles are on positions (sites) of a finite subset V C Z d , and any particle has 
either spin +1 or —1. So let A = { — 1, +1} and let E = A v . Any element of E shall be called 
configuration of the system. Configurations will be denoted by letters rj, £, (. Any element i £ V 
is called a site. For any site i let rf be the modified configuration r] l {i) = —T}{i) : rf{j) = r](j) for 
all j 7^ i. We associate to any site i and any configuration r\ a spin-flip-rate Ci{rf) > 0. Roughly 
speaking, site i will change its spin at rate Ci(t]) whenever the overall configuration of particles 
is rj. We suppose that 

sup Ci(ri) < Mi 
v 

for some constant Mj. Then the associated interacting particle system is a Markov process X 
on E having generator 

i/(i) = E^)i/w-/wi- 

Let Vi{k) = {j G 7L d \ < \\j — i\\ < k}, where ||j|| = Yut=i \^ u \ ^ s ^ e usua l -^l-norm of Z d . In 
Galves et al. [11] the following criterion for recurrence of the system has been given : 

Proposition 6.1 (Theorem 3 of Galves et al.) There exists a sequence Xi(k),k > —1, as- 
sociated to the spin-flip rates q such that : If 

supV \Vi(k)\Xi{k) < 1, (6.50) 

then the process is recurrent. If the above sum is strictly less than 1, then the process is uniformly 
exponentially ergodic. Here, \Vi{k)\ is the number of sites belonging to Vi(k). 

Remark 6.2 The above theorem relies on the construction of a backward dual process 
0, such that for any site i G V, CT 1 denotes the set of sites at time — s that have to be known 
in order to determine the spin of site i at time 0. The cardinal of can be compared to a 
classical branching process in continuous time, having reproduction mean ^2 k>0 |V^(A;)|Aj(A;) < 1, 
and thus being subcritical. We refer the reader to Galves et al. (\TT$) for the details. 

The sequence \{k) can be constructed explicitly, we refer the reader to Galves et al. (2008) 
for the details. In [llj . we were mainly concerned with the issue of perfect simulation of the 
invariant measure /1 of the process, under condition (16.501) . based on the precise knowledge of 
the spin-flip rates q. Suppose now, that we are in the converse situation, observing the process 
over some time interval [0,t], without knowledge of the spin-flip rates, and that we want to 
deduce information about the associated invariant measure (and for example on the associated 
spin-flip rates, a posteriori). This is the classical situation of statistical inference. 
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Note that assumption 2.1 is satisfied and that the process is strongly Feller. 

We apply Theorem 14.21 with rj — \. Since the process is uniformly exponentially ergodic, any 
constant function is special. In particular, we have for any bounded function /, since / — //(/) 
is special, that 

p(\~J*f(X s )ds-Li(f)\>x S j <ce- ct ^ + R t (x), (6.51) 

for some constants c, C. This is a first step on the way towards estimating spin-flip rates and 
associated interaction schemes for interacting particles. 

7 Final remarks 

Our results obtained in the positive recurrent case (Theorem I3.1[ item 1.) have to be compared 
to the results obtained by Loukianova et al. (|18j). The results stated in the present paper 
are different from the situation studied in Loukianova et al. (2009) since we consider centered 
functions that are special, i.e. satisfying sup^E^, J^ 1 \f\(X s )ds < oo. Additive functionals built 
of special functions admit exponential moments independently of the degree of recurrence of 
the underlying process, and this is a very special feature about special functions that does in 
general not hold for other functions. 

Namely, if / is any positive bounded function having compact support, then / is certainly 
special, but f '■— f — /•*(/) i* 1 general won't (constants are special only if the underlying process 
is uniformly ergodic). 

Note however that our method works for any function / such that exponential moments ( I4.27P 
are finite. Thus if, for A sufficiently small, 

E n (e^ R2 - Rl) ) + E w (e Ai?1 ) < oo, (7.52) 

then our result holds also for any function f = f — fi(f), where / is positive, bounded and of 
compact support. It will be the subject of some futur work to give necessary conditions for 
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